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; SPATIAL PAFPERN AND INTERFERENCE 
IN PINON-JUNIPER WOODLANDS OF NORTHWEST COLORADO 


Charles W. Welden!*, William L. Slauson!, and Richard T. Ward! 


ABSTRACT.—The local spatial arrangement of the coniferous trees Pinus edulis and Juniperus osteosperma was 
mapped in two woodland stands and measured in two shrub-dominated stands in the semiarid Piceance Basin of 
northwest Colorado. In the woodlands, small trees were olten clumped, while medium and large trees were cither 
randomly or uniformly dispersed. Significant regressions were obtained between a tree’s basal area or Canopy area and 
the area ofits Dirichlet domain (the region closer to it than to any other tree). Both findings from the woodland stands 
accord with results obtained by other workers in other vegetation. Like earlier workers, we interpret these patterns to 
indicate density-dependent mortality and densitv-dependent depression of growth rates among the trees in the 
woodlands. In contrast, the trees in the shrub-dominated stands are located at random with respect to each other. 
Hlowever, they are strongly associated with shrub cover. Apparently, tree seeds arrive in these stands primarily by 
long-distance dispersal, and the establishment of seedlings is more likely in the shade of shrubs. 


Since plants are sessile and their growth is 
plastic, their arrangement in space and their 
sizes can reflect the history of their interac- 
tions with each other and with the environ- 
ment. With long-lived, slow-growing plants, 
studying pattern may be the only feasible 
way to discover which processes and inter- 
actions are important in determining commu- 
nity structure. 

We used some of the methods compared 
by Goodall and West (1979) to study the local 
spatial arrangement (pattern) of the small co- 
niferous trees Pinus edulis and Juniperus 
osteosperma in four stands in the semiarid 
Piceance Basin of northwest Colorado. Our 
goals were twofold. First, we wished to de- 
termine whether the differences between 
methods Goodall and West (1979) detected in 
artificial populations are borne out in more 
complex real populations. Second, we wished 
to infer the processes that influence the estab- 
lishment of seedlings and the growth and mor- 
tality of plants. 


STUDY AREA 


The Piceance Basin occupies about 3000 


ts 2 - . . 
km* in Garfield and Rio Blanco counties of 


northwest Colorado. Elevations range from 
1707 to 2743 m (Tiedeman and Terwilliger 
1978). The climate is semiarid with average 


annual precipitation ranging from 28 em in the 
northwest to 63.5 cm in the southeast. About 
half of the annual total falls as snow and most 
of the remainder as rain in late-summer thun- 
derstorms. In the short term, precipitation is 
unpredictable and variable (Wymore 1974). 

The average annual temperature is 7 C at 
1825 m (the elevation of the only permanent 
weather station in the basin), with a minimum 
monthly average in January of —5.9 C and a 
maximum monthly average in July of 20.3 C. 
The average annual temperature decreases by 
approximately 0.85 C for every 100 m increase 
in elevation. Both temperature and precipita- 
tion are strongly influenced by local topogra- 
phy (Wymore 1974). 

We studied the spatial patterns of Pinus 
edulis Engel. and Juniperus osteosperma 
(Torr.) Little (pifon and Utah juniper). 
Nomenclature follows Goodrich and Neese 
(1986). P. edulis and J. osteosperma are small 
coniferous trees common throughout the 
western United States, where they form 
mixed stands, often with an understory of 
scattered grasses, forbs, and shrubs. They 
commonly attain heights of 6-8 m, and both 
reproduce by seed. P. edulis usually possesses 
a single stem, while J. osteosperma is often 
multistemmed. 

The vegetation of the basin includes 
shrublands and woodlands of various floristic 
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compositions. Pinon-juniper woodlands (as 
deseribed in Tiedeman and Terwilliger 1978) 
have open eanopies dominated by P. edulis 
and J. osteosperma and occur on broad, flat 
ridge tops at elevations between 1890 m and 
2170 m, where soils are shallow, roeky, light 
brown, sandy loams (Entisols). Shrublands 
dominated by Artemisia tridentata Nutt. 
(sagebrush flats) often occur on the same 
ridges as do pinon-juniper woodlands, at 
roughly the same elevations, but where soils 
are finer and deeper. Where pinon-juniper 
woodlands abut sagebrush flats, zones of in- 
termediate vegetation are often found. In 
these intermediate areas, the vegetation is 
dominated by Artemisia, with small, scat- 
tered individuals of P. edulis and J]. osteo- 
sperma. Few of the trees overtop the shrubs. 

We studied two pinon-juniper woodlands 
(stands A and B), which were dominated by 
mature P. edulis and J. osteosperma, with 
little shrub understory. The canopies in these 
stands are not closed, but individual cano- 
pies sometimes abut or overlap. It is known 
from others (Fowells 1965) and from personal 
observation that the roots of these trees usu- 
ally extend beyond the canopy. Thus, neigh- 
boring trees which do not seem to be compet- 
ing for light may nonetheless be competing 
belowground for water or nutrients. These 
stands lie at elevations of 2164 m and 1890 m, 
which approximate the elevational limits of 
this vegetation in the basin. Stand A slopes 
1.5° and faces to the northwest (N62°W). 
Stand B slopes 3.0°, facing to the north-north- 
west (N22°W). 

Stands C and D are intermediate between 
pinon-juniper woodlands and sagebrush flats. 
None of the trees in these stands is as large as 
the largest trees in the pifon-juniper wood- 
lands, although many bear cones and are thus 
sexually mature. These stands oceupy ridge 
tops at elevations of 2164 mand 1981 m. Stand 
C slopes 4.5°, facing west (NSO°W), and stand 
D slopes 6.5°, facing north (N5°W),. 


METHODS 


Goodall and West (1979) reviewed pattern 
methods based on analyses of artificial popula- 
tions. They compared the statistical powers 
of the methods, that is, the probabilities of 
rejecting a false null hypothesis. With large 
samples, all the tested methods gave results 


COW. WELDEN ETAL. 


[Volume 50 


reflecting the true dispersion pattern of artifi- 
cial populations, with powers approaching 
100%. With smaller samples, however, 
methods differed in power. We used those 
having the greatest power with small samples: 
the variance/mean ratio (Clapham 1936) 
among quadrat methods, and the indices of 
Hopkins (1954) and Pielou (1959, 1960, 1961) 
among distance methods (see descriptions be- 
low). We also compared the frequencies of 
quadrats containing exactly 0, 1, 2, ... plants 
with the expected Poisson distribution by a 
chi-squared goodness-of-fit test. 

In addition to these methods, we included 
a measure of pattern that uses information 
not only about the locations of plants but also 
about their sizes. The Dirichlet domain (or 
Thiessen or Voronoi polygon) of a plant com- 
prises all the points closer to that plant than to 
anv other (Honda 1978, Jack 1967, Mead 
1971, Mithen, Harper, and Weiner 1984), Its 
size thus represents the area more easily ae- 
cessible to the plant than to its neighbors and 
mav represent the amount of resources cap- 
tured or sequestered by a plant, or potentially 
more available to it than to its neighbors. This 
in turn may influence the plant's growth and 
fitness and indicate what effect, if any, its 
neighbors have on it. To detect whether this is 
the case, we regressed the areas of plants’ 
Dirichlet domains on the sizes of the plants. 

The variance/mean ratio test (Clapham 
1936) is based on the expectation that, in a 
randomly dispersed population, the fre- 
quency distribution of quadrats containing ex- 
actly 0, 1, 2, 3,... individuals approximates 
the Poisson distribution. One property of this 
distribution is that its mean and variance are 
equal, and their ratio therefore unitv. The 
distribution of this ratio in large samples is 
approximately normal, with a mean of | anda 
standard deviation of (2/n—1)"” (Goodall and 
West 1979), where » is the sample size (num- 
ber of quadrats). In regularly dispersed popu- 
lations the ratio is less than 1, in aggregated 
ONeCSs 21 Cater. 

Hopkins’ (1954) index A is based on the 
expectation that, in a randomly dispersed 
population, the average distance from ran- 
domly located points to the nearest plant 
equals the average distance between plants 
and their nearest neighbors. Hopkins pro- 
posed the ratio of these two averages as his 
index: 
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where P, and I, are the sums of equal numbers 
of distances from random points to the nearest 
plant and from randomly selected plants. to 
their nearest neighbors, respectively. In a 
randomly dispersed population, the expected 
value of A is 1, and for large samples its fre- 
quency distribution is approximately normal. 
Values of A larger than | indicate aggregation, 
less than | regularity. 

Pielou (1959, 1960, 1961) developed two 
distance methods to measure pattern. The 
first uses a sample of distances from randomly 
located points to the nearest plant and an 
independent estimate of plant density. From 
these a statistic, alpha,, can be calculated as 
follows: 


py? 


alpha, — pi(D)omega, 


where D is the density of the plants, omega, 
is the mean squared point-to-plant distance, 
and pi is the trigonometric constant. 

The second method (Pielou 1960) uses a 
sample of distances from randomly chosen 
plants to their nearest neighbors. A statistic 
alpha, is calculated in the same way as alpha, 
substituting the mean squared plant-to-plant 
distance for the mean squared point-to-point 
distance. Pielou (1959) provides tables of con- 
fidence intervals and significance levels for 
values of alpha, and shows how they may be 
used to interpret alpha, (Pielou 1960). 

We mapped the location of each Pinus 
edulis and Juniperus osteosperma 10cm tall or 
taller in parts of stands A and B. The mapped 
area in stand A was 2250 m: in stand B it was 
2500 m>. We checked the accuracy of the maps 
by comparing plant-to-plant distances calcu- 
lated from map coordinates to the same dis- 
tances measured in the field. The greatest 
difference was about 10 cm. 

We classified plants into three height- 
classes. Small plants were 10 cm to I m tall, 
medium plants between 1m and 3 m tall, and 
tall plants were taller than 3 m. The tallest 
trees in our stands were about 5 m tall. Small 
plants were not mapped in about one-third 
of stand A. 

For each P. edulis in these stands we mea- 
sured one canopy diameter in an arbitrary 
direction and estimated the area of its canopy 
as if it were circular. The living canopies of 
J. osteosperma were often interrupted by 
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dead branches. We measured the living 
portions of their canopies and summed the 
areas estimated from these. Basal areas were 
calculated for both species from stem diame- 
ters measured at ground level. For multi- 
stemmed plants, the basal areas of all living 
stems were summed. 

We measured the dispersion patterns of the 
plants on these maps, using both quadrat and 
distance methods. Small plants were sampled 
with quadrats 2.5 m on a side (in map scale), 
medium and tall plants with quadrats 5 mona 
side. Quadrats were placed at the intersec- 
tions of a regular grid of lines 5 scale-meters 
apart; thus every point on the map was in- 
cluded in exactly one quadrat of a given size. 
There were 100 large and 400 small quadrats 
in stand B. Stand A was more irregular, en- 
compassing 90 large and 230 small quadrats. 

The spatial dispersions of each size class 
and species were measured separately and 
pooled. That is, the null hypothesis of random 
spatial dispersion was tested by five indices 
for small P. edulis, small J. osteosperma, all 
small plants, medium P. edulis, medium J. 
osteosperma, all medium plants, tall P. edulis, 
tall J. osteosperma, all tall plants, medium and 
tall P. edulis combined, medium and _ tall 
J. osteosperma combined, and all medium 
and tall plants combined. 

We constructed Dirichlet domains (Honda 
1978, Jack 1967, Mead 1971, Mithen, Harper, 
and Weiner 1984) for the plants by drawing 
lines connecting each plant to its immediate 
neighbors, and then constructing perpendic- 
ular bisectors of these lines (Fig. 1). Note that 
we did not weight the distance from a plant to 
the bisector by the size of the plant, and thus 
there is no necessary correlation between the 
size of a plant and the size of its Dirichlet 
domain. We estimated the areas of the Dirich- 
let domains by cutting the polygons from the 
maps and weighing them. We regressed the 
areas of the Dirichlet domains on the basal 
areas, and separately on the canopy areas, of 
their plants. Regressions on basal areas were 
compared to regressions on canopy areas, 
with and without logarithmic transformation, 
by graphical analysis of residuals. 

In stands C and D we located every P. edulis 
and J. osteosperma 10 cm or more in height 
within a square 50 m on a side, noting whether it 
had become established under a plant canopy or 
in the open, based on observations of each 
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Fig. 1. Construction of Dirichlet domains. 





a. Draw line segments connecting a focal plant to its 
neighbors. 

b. Draw the perpendicular bisector ofeach line segment. 

c¢. The Dirichlet domain is the region closer to the focal 
plant than all the perpendicular bisectors. 

d. Repeat for each plant. The Dirichlet domain of each 
plant is the region closer to it than to any other plant. 


tree’s association with living or dead shrubs. 
We measured total plant cover of all species 
with two 50-m line intercepts. The association 
of P. edulis and J. osteosperma with plant 
cover was tested by a chi-squared test. We did 
not map these stands but measured distances 
between neighboring trees in the field. We 
used Pielou’s alpha, to describe the spatial 
dispersion of the two tree species. 
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RESULTS 


Table 1 shows the number of P. edulis and 


J. osteosperma in each stand and the corre- 


sponding numbers per hectare. Table 2 shows 
the five dispersion indices for the trees in 
stands A and B, and Table 3 the interpreta- 
tions of these values. In the woodland stands 
(A and B) small plants tend to be clumped, and 
larger plants tend to be randomly or uniformly 
dispersed. The sequence from clumped to 
random to uniform is violated in only three 
instances (asterisks in Table 3). These viola- 
tions may be the result of chance, since the 
tests for significance were all set at the 5% 
level and some spurious results are expected 
among such a large number of separate tests. 
All log-log transformed regressions of 
Dirichlet domain areas on plant canopy 
areas and basal areas in stands A and B are 
significant at the 5% level, except for that of 


J. osteosperma in stand A (Table 4, Fig. 2). 


These regressions show that, on average, 
larger plants have larger Dirichlet domains 
and are correspondingly farther from their 
neighbors. The Dirichlet domains of small 
plants are more variable in area than those of 
larger plants. Logarithmic transformation of 
both variates improves the distribution of 
variates and residuals and produces reason- 
able conformity with the assumptions of re- 
gression, but it does not change the signifi- 
cance ol the regressions. These results are 
similar to those of regressing the distance be- 
tween a pair of neighboring plants on the sum 
of their sizes (Welden 198-4, Welden, Slausen, 
and Ward 1988, cf. Fuentes and Gutierrez 
19S1, Gutierrez and Fuentes 1979, Nobel 1981, 
Phillips and MacMahon 1981, Pielou 1960, 


TABLE |. Stand censuses, divided by height categorics (10cm = small < 1m = medium < 3m & tall) and by species. 


In parentheses are numbers per hectare. 








Small Medium tall Total 
plants plants plants 
Stand A Pinus edulis 61 (432) 36 (229) 67 (274) 184 (753) 
Juniperus osteosperma 26 (15-4) [2 (49) 39 (160) 77 (315) 
Stand BP. edulis 5S (352) 32 (128) 1) (44 131 (529) 
J. osteospermia 86 (3-E4) 9 (36) ‘EL (16-4) 136 (5-14) 
Stand CP. edulis 56 (224) 22 (SS) Lets 79 (316) 
J. osteosperma a 2) 7 12s) 0 14 (56) 
Stand DP. edulis 17 (LSS) ou) (120) 4 (16) SI (324) 
J. osteosperma SLOr36) 46 (1S-f) [7 (6S) 97 (38S) 
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TABLE 2. Values of dispersion indices in stands A and B. fndices are identified in the text and these values are 
interpreted iu Table -f. A dash indicates that the test could not be performed. 


Distance methods 


Quadrat methods 





Stand alpha, alpha, A chi? var/mican 
A B A B A BB A B AN B 

P. edulis 

Simall 0.93 ye! 2 0.S4 0.83 e271 eee 25.0 ear [.9S 

\ledium 0.98 123 L.15 0566 0.85 leo3 SOF _ lai 1.60 

Medium and tall  b.00 9 0.74 1.04 0.97 0.96 1.38 2.90 2.16 [4S 1.28 

Tall O.S4 0.358 O.SS 0.9] 0.95 Ol Qo — 1.10 0.S9 
|. osteospernia 

Small 1.03 1.36 1.04 1207 0.99 L237 — 29.4 Lea is: 

Medimu 0.6] 1.06 0.54 0.64 O72 iets = 1.04 0.86 

Medium and tall 0.73 0.93 1.09 1.29 0.78 0.73 ( ||P 0.96 0.86 

Tall i arcs eS: E.1S 139 Ov7a 0.54 1.10 0.97 0.82 
Species combined 

Small L225 1.29 0.79 0.S6 1.58 Heat) 1S. ] S.1S 1.62 1.65 

Medium ay 1.58 0.96 0.54 eas 2O4 225 10:9 ee 22) 

Medium and tall 1.06 0.86 (). 94 L.14 ele 4) 92 top O91 bot 1.19 

Tall ().92 0.74 0.98 Tals 0.78 0.59 0.60 ahs 1.08 0.95 


TABLE 3. Pattern analyses of stands A and B. C indicates that the plants arc clumped, R that they are randomly 
dispersed, U that they are uniformly dispersed. AH indicated nonrandom dispersions are significant at the 5% level. A 
dash indicates that the test could not be performed. Asterisks denote contradictions to the general treud of C — R— U 


with increasing plant size. 





Distance methods 


Stand alpha, alpha, 
A B A B 

P. edulis 

Small R R R R* 

Medium R R R © 

Medium and tall R U R R 

Tall R U R R 
J. osteosperma 

Small R C R R 

Medium La) R R R 

Medium and tal] U R R R 

Tall U U R U 
Species combined 

Small Cc Cc R RK 

Medium R Cc R CG 

Medium and tall R R R R 

Tall R U R R 


Quadrat methods 


A chi? var/mean 
A B A B A 3 
R R* C C C C 
R ce C — ( fc 
R R R R ( R 
R R KR — R R 
R R — CG C C 
R R — — R C 
R R R R R R 
R U R R R R 
R (C (© C C c 
R ee hh GC G, C 
R RK R R c R 
R U RK R R R 





1961, Yeaton and Cody 1976, and Yeaton, 
Travis, and Gilinsky 1977). 


Plant cover (primarily of Artemisia ) in stand 


C was approximately 20%, and about 96% of 


the P. edulis and about 71% of the J. osteo- 
sperma had become established under plant 
canopy. Plant cover in stand D was about 
18%, and about 93% of the P. edulis and about 
87% of the J. osteosperma had become estab- 
lished under plant canopy. The probability 
that establishment of P. edulis or ]. osteo- 


sperma is random with respect to plant cover 
is less than .0O1 in every case. The pattern 
statistic alpha, (Pielou 1960) showed no sig- 
nificant deviations from random dispersion 
among P. edulis or J. osteosperma in stands 


Cand D. 
DISCUSSION 


Pielou (1959) and Goodall and West (1979) 
show that distance methods are more sensi- 
tive to uniformity and quadrat methods are 
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TABLE 4. Coefficients of log-log transformed regressions of Dirichlet domain area on canopy and hasal areas. 
Significance is the probability of such data ifthe true slope and r equal zero. 


Species Stand n 


P. edulis 
Independent variable 
In (canopy area) 
In (basal area) 
In (canopy area) B 33 
In (basal area) 


J. osteosperma 


In (canopy area) A 27 
In (basal area) 
In (canopy area) B 31 
ln (basal area) 
Species combined 
ln (canopy area) A 125 
ln (basal area) 
In (canopy area) B 6-4 


In (basal area) 


15 


@ 


In(Dirichlet Area) (cm*) 





15 


o 10 > 
In(Canopy Area) (cm ) 


Fig. 2. Regression of Dirichlet domain area on canopy 
area of pifons in stand A. Both variates have been trans- 
formed to their natural logarithms. 


more sensitive to clumping. This is borne out 
by Table 4, where it can be seen that the 
quadrat methods never detected uniform dis- 
persion while the distance methods did. The 
distance methods, on the other hand, failed to 
detect clumping in several cases where it was 
detected by the quadrat methods. 

The trees in the woodland stands (A and 
B) appear to be interfering (sensu Harper 
1961, 977) with one another, either by com- 
petition or by allelopathy. The trend from 
clumped to random to uniform dispersion 
with increasing plant size suggests densitv- 
dependent mortality. Density-independent 
mortality ina clumped population might con- 


93 


ae Y-intercept Slope Significance 
0,056 10.12 0.13 0.019 
0.052 11.05 0.11 0.024 

on 9.15 0.26 0.002 
0.191 6.5] 0.21 0.011 
0.031 12.50 —(.06 0.377 
0.020 12.04 —0.04 0.479 
0,367 3.4] 0.23 0,000 
0.352 6.54 0.18 0.000 
0.039 10.49 0.10 0.027 
0.047 }1.15 0.09 0.015 
0.333 Ooo 0,25 0.000 
0.265 6.58 0.18 0.000 


ceivably reduce sample sizes in successively 
larger size-classes until the clumping is no 
longer detectably different from a random dis- 
persion, but it seems unlikely that it could 
produee a uniform dispersion (Phillips and 
MacMahon 1981). 

The significant regressions of Dirichlet 
domain area on plant size indicate density- 
dependent mortality or density-dependent 
suppression of growth, or both. We envision 
two processes leading to this result. First, 
plants that become established farther from 
preexisting neighbors become larger because 
they have access to more unexploited (or un- 
sequestered) resources. Second, established 
plants prevent the establishment of new 
neighbors nearby, or impede their growth, 
because they have exploited (or sequestered) 
most of the resources in their neighborhoods. 

Mithen, Harper, and Weiner (1984) found 
significant positive relationships between 
Dirichlet domain area and plant dry weight in 
even-aged greenhouse populations of Lap- 
sana communis L. Although the conditions of 
their experiments are different (particularly 
since their plants germinated svnehronously), 
their interpretations of their results are simi- 
lar to ours here. 

Pielous (1959, 1960) method did not de- 
tect any deviation from random spatial ar- 
rangement in stands C and D. Ekowever, 
both tree species are significantly associated 
with plant cover. We presume that these trees 
became established after long-distance dis- 
persal (> 100 im) from nearby woodlands. The 
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significant interaction in these stands is evi- 
dently not interference between neighboring 
trees, but amelioration of abiotic stress under 
the eanopies of preexisting plants. Fowells 
(1965) reports that P. edulis requires shade 
early in its development. 

Our evidence for these interpretations is 
circumstantial. However, given the long lives 
and slow growth of these plants, and the vary- 
ing physical environment of the study area, 
sueh evidence may be the most informative. 
These pattern methods integrate the effects of 
environment and biotic interactions over the 
life spans of the plants, a time seale not usually 
accessible to more mechanistic methods. 

All our inferences of processes leading to 
the present pattern require further examina- 
tion. Although J. osteosperma has been re- 
ported to produce allelochemicals (Jameson 
1971), experiments should be done to deter- 
mine whether allelopathie effects occur under 
the conditions and in the soils of the Piceance 
Basin, and more field studies are needed to 
determine whether establishment occurs 
more often near neighbors or far from them. 
The dynamic behavior of the various pattern 
indices and regressions should be explored 
under conditions of densitv-dependent and 
density-independent mortality. 
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